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Abstract. The problem of counting monomer-dimer coverings of a finite patch in a lattice 
is a longstanding problem in statistical mechanics. It has only been exactly solved for the 
special case of dimer coverings in two dimensions (U], IHl)- In earlier work, Stanley (|H]) 
proved a reciprocity principle governing the number N{m, n) of dimer coverings (also known 
as perfect matchings) of an m by n rectangular grid, where m is fixed and n is allowed to 
vary. As reinterpreted by Propp (0), Stanley's result concerns the unique way of extending 
N{m, n) to n < so that the resulting bi-infinite sequence, Nim, n) for n £ Z, satisfies a linear 
recurrence relation with constant coefficients. In particular, Stanley shows that N{m, n) is 
always an integer satisfying the relation N{m, —2 — n) = tm,nN(m, n) where em,n = 1 unless 
TO = 2 (mod 4) and n is odd, in which case tm,n = —1. Furthermore, Propp's method 
was applicable to higher-dimensional cases, such as the dimer model on an mi-by-m2-by-n 
grid. This paper discusses similar investigations of the number M(m,?i), of monomer-dimer 
coverings, or equivalently, (not necessarily perfect) matchings of an m by n rectangular grid. 
We show that for each fixed m there is a unique way of extending M{m, n) to n < so that 
the resulting bi-infinite sequence, M(m, n) for n S Z, satisfies a linear recurrence relation with 
constant coefficients. We show that M(m, n), a priori a rational number, is always an integer, 
using a generalization of the combinatorial model offered by Propp. Lastly, we give a new 
statement of reciprocity in terms of multivariate generating functions from which Stanley's 
result follows. 



1. Introduction 

1.1. Background. The problem of counting the monomer-dimer coverings of a finite patch of 
a lattice has been examined for many years in the field of statistical mechanics (see Q), and has 
applications in biology, chemistry and physics. The closely related dimer problem, where the 
number of monomers is zero, was exactly solved for two dimensional lattices by Kasteleyn ([3]) 
and Temper ley and Fisher ([HI)- While we cannot offer an exact solution to the monomer-dimer 
problem, even in two dimensions, we can describe a symmetry property that must be satisfied 
by exact solutions to various restrictions of this problem, in which one considers a sequence of 
subgraphs of a lattice, where the successive subgraphs grow in one direction while their size in 
all other directions remains fixed. 

In particular, we consider families of subgraphs of lattices indexed by the natural numbers 
(n = 1, 2, 3, . . .), and count the number a„ of monomer-dimer coverings of the nth subgraph. The 
numbers a„ satisfy a linear recurrence with constant coefhcients and therefore are given by an 
exact formula of polynomial-exponential type; one can thus in each individual case write down 
an exact formula for a„, at least in principle, and in particular, one can study symmetries of the 
resulting function of n when n is no longer restricted to being a natural number. Alternatively, 
when n is a negative integer, one can often find a„ by working directly with the recurrence 
relation satisfied by the sequence. Surprisingly, we find that the new numbers are related, 
term by term, to the original sequence of numbers; for instance, in some cases an ~ a-n-2- 
Richard Stanley's book 7 , in the context of rational generating functions, devotes an entire 
section to exploring the relationships (called reciprocity relationships) between positively- and 
nonpositively-indexed terms of a sequence. The nonpositively-indexed terms may even have a 
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nice combinatorial meaning on their own, as seen in the case of Ehrhart reciprocity, for example. 
In this paper we describe a relationship between terms of certain integer sequences {a„}, n € Z, 
as well as give a combinatorial interpretation to the "un-natural" terms. 

1.2. The Problem. Rather than use the language of statistical mechanics (monomers and 
dimers on lattices), we will refer to matchings of graphs where we define a matching of a graph 
G = {V, E) to be a collection of edges of G, no two of which share a vertex, together with all 
the vertices of G that are not incident with those edges. ^ Edges correspond to dimers, isolated 
vertices correspond to monomers. A matching is called perfect if it is composed entirely of edges, 
as in the dimer problem. If we consider the number of matchings of a rectangular grid-graph of 
fixed height m, and varying width n, we obtain an integer sequence, {a„}. For example, when 
m ~ 1, we obtain the Fibonacci sequence: 

1,2,3,5,8,... 

This sequence satisfies a linear recurrence, namely a„ — an~i + an-2- So we can always obtain 
the value of a term based on the two terms to the left of it. But similarly, we can obtain the 
value of a term by the two terms to the right of it, e.g. 3 = 8 — 5, or a„ = a„+2 — dn+i- In this 
way we can extend the Fibonacci sequence to the left and get values for a„ when n < 0: 

... - 8, 5, -3, 2, -1, 1, 0, 1, 1, 2, 3, 5, 8, . . . 

So we see that we now have a doubly infinite integer sequence, and that it is symmetric up to 
sign, i.e. a„ = ±a_„_2. Some natural questions arise: Does this symmetry mean something? 
Do the values of a„ for n < 0, being integers, count something (at least up to sign)? If so, can 
we extend the result to larger values of ml 

The answer to all the questions is yes. In 5 , Propp considered integer sequences N(rn,n), 
generated by perfect matchings of rectangular grid graphs of fixed height, and how they extend 
to values for n < 0.^ He came up with a unified combinatorial object and a way of counting 
signed matchings of these objects that allowed him to show that any half-infinite sequence given 
by the number of perfect matchings of a fixed-height grid graph extends to a bi-infinite sequence 
that does three things: 

• It satisfies a linear recurrence relation of finite degree with constant coefficients. 

• The numbers obtained when going backwards are unique and are always integers. 

• The bi-infinite sequences have a special kind of symmetry, or reciprocity, stated roughly 
as \N{m, n)\ — \N{m, ~n — 2)| for any to > 0. 

Matchings in general (of the graphs we will classify) satisfy the same sort of linear recurrence 
(as we will explain in section 3). But even given that a covering problem satisfies a linear 
recurrence, integrality is not ensured when running the recurrence backwards.* For a simple 
example, consider covering a 2-by-n grid with monomers and dimers, but only allow the dimers 
to be vertical. (It is easy to see that the number of coverings is 2", but in keeping with the 
spirit of this article, we ignore the exact formula and work with a linear recurrence instead.) 
The number of such coverings, a„, is governed by a„ = 2a„_i. Upon observing that ai = 2, 
we can generate the half- infinite sequence 2,4,8,16, To work the sequence backwards, we 



Ehrhart reciprocity describes a relationship between the number of lattice points found in a closed rational 
polytope, and the number of lattice points found in its interior. See 0. 

^Another equivalent terminology would be to speak of tilings of planar regions where the tile set consists of a 
domino (2 by 1 rectangle) and a monomino (1 by 1 rectangle). 

■^His first nontrivial case, N{2, n) is also the Fibonacci sequence. 

priori the numbers obtained from running a linear recurrence in reverse need only be rational. If the 
recurrence is of the form a„ = cia„—\ + cgfln— 2 + ■ ■ ■ + Cfea„_j., then to push the sequence backwards we solve 
tor the term oi smallest index to get a„-fc = 
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reverse the recurrence by writing it as a. 



^a„+i and then generate the bi- infinite sequence 




For n = —1 we cease to observe integrahty. However in the case of monomer- dimer tihngs, 
we can guarantee integrahty of the sequences by showing they can be generated by counting 
matchings of certain graphs for all values of n. 

The objects we will present not only ensure integrality, but also give a direct proof of the 
paper's main idea: a reciprocity statement for general matchings. This claim would appear to 
be the least obvious judging by the numerical evidence. Indeed, consider the integer sequence 
M(2, n), the number of matchings of a 2-by-n rectangular grid graph. The sequence is governed 
by the recurrence 



with initial 'natural' values M(2, 1) — 2, M{2, 2) = 7, and M{2, 3) = 22, allowing us to generate 
the bi-infinite sequence 



Although we get integers when going backwards, we seem to lack the symmetry that Stanley 
and Propp observed for perfect matchings. But a premonition of symmetry can be seen if one 
considers the parity (odd versus even) of each term; one finds that M(2, n) and M(2, — 2 — n) 
are congruent mod 2. 

Why there should be symmetry (or symmetry mod 2) in integer sequences as above is ex- 
plained by the answer to the question: what do the numbers, M(m, n) for n < 0, count? Propp 
had an answer for the case of perfect matchings, but it does not directly translate to general 
matchings. Section 2 outlines the dilemma and gives a solution by way of objects called "empty 
vertices" . In weights were not given to isolated vertices, simply because there were none. But 
when examining matchings as we have defined them (a monomer-dimer model) , isolated vertices 
do have weights, giving rise to the need for empty vertices to explain the observed phenomenon 
and tell the "proper" combinatorial story. By incorporating empty vertices we will obtain a 
combinatorial model whose (signed) number of signed matchings is M{m,n) for all n. Looking 
at the signed matchings we will be able to see the symmetry hidden in the terms above: 

... 121 - 107, 41 - 30, 12 - 10, 5 - 2, 1 - 1, 1, 0, 1, 1 -f 1, 5 + 2, 12 + 10, 41 + 30, 121 + 107, . . . 

This also explains the symmetry mod 2. 

The properties exhibited by the numbers M(m, n) and the objects they count can also be 
seen in terms of the generating function Fm{t, x, y, z) = Yl^=i fnix, y, z)t". Here x is the weight 
given to horizontal edges, y is the weight given to vertical edges, and z is the weight given to 
vertices. The polynomial fn{x, y, z) encodes all of the matchings of an to by n grid graph. We 
will show 



The second item in the list shows how perfect matchings may be extracted from general match- 
ings. The last item in the list is a heuristic expression of reciprocity that will be made precise 
in section 4.2. 

The theory of empty vertices is applicable to more general situations than simply matchings of 
rectangular grid graphs. We will prove results for matchings of graphs of the form G x P„, where 
G is an arbitrary finite graph and P„ is the path graph of length n. This is the "box-product" 
of the two graphs: a vertex of G x P„ is a pair (u, v) consisting of a vertex of G and a vertex of 
P„, and two such vertices (m, w), (u',w') are connected by an edge iff either u ~ u' and v^v' are 
adjacent in P„ or v — v' and u, u! are adjacent in G. Wc will refer to such graphs as "generalized 



M(2, n) = 3M(2, n - 1) + A/(2, n - 2) - Af (2, n - 3) 



. . . 14, 11, 2, 3, 0, 1, 0, 1, 2, 7, 22, 71, 228, . . . 



(1) P„(i,l,l,l) = E^^iM(TO,n)i" 

(2) P„(i,l,l,0) = Er=i^K")*" 

(3) Fm{t,x,y, z) ~ ~F„i{l/t,x, -y, -z) 



4 



NICK ANZALONE, JOHN BALDWIN, ILYA BRONSHTEIN, AND T. KYLE PETERSEN 



rectangles". Notice that an to by n rectangular grid graph is G x P„ where G — Pm- If M{n) 
is the number of matchings of G x P„, we can show that there are objects (signed graphs) for 
which the (signed) number of matchings is M{—n), i.e., we can give a definition for G x 
Reciprocity asserts that the (signed) number of signed matchings of G x P_2-n is equal to (or 
equal to the negative of) the (unsigned) number of matchings of G x P„. Throughout most of 
the paper, our proofs will be for generalized rectangles, but examples will usually involve grid 
graphs. 

2. Signed Graphs and Signed Matchings 

Let us try to deduce the nature of to by n rectangular grid graphs for all n, not concerning 
ourselves with matchings at all. In Figure^ denote the 2 by 3 grid graph by G(2, 3). The graph 
G(2, 4) is obtained from G(2, 3) by adding two horizontal edges, a vertical edge, and two vertices 
as shown. Likewise, we can obtain G(2, 2) by removing the same set of edges and vertices from 
the right of G(2, 3). Then we get G(2, 1) and G(2, 0) inductively. But what are we to make of 
G(2, 0)? It has no vertices, and two "anti-" horizontal edges. How many matchings should this 
have? What about G(2,-l)? G(2,-2)? 



G(2,3) 



G(2,3) 



G(2,2) 



G(2,l) 



G(2,4) 



G(2,2) 



G(2,l) 



G(2,0) 



G(2,0) 



■o- 



■O- 



G(2,-l) 



Figure 1. Moving from G(2,3) to G(2,4) or G(2,2) 

The definitions that Propp created for G(to, n) with n < look quite similar to the pictures 
in Figure n except that his horizontal lines are solid and there are vertices on the endpoints of 
the horizontal edges where there are none in Figure ^ (see Figure . 



Figure 2. Propp's model for G(2, -1) 



Notice that in Figure ^ the graphs for G(2, 0) and G(2, —1) have some horizontal edges that 
are incident with only one vertex. The meaning of these missing vertices, these places that are 
"empty" of vertices, is described below. Let G — (V, E) be a graph in the usual sense, i.e. a set 
of vertices and edges between them, except that there is more than one kind of vertex and more 
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than one kind of edge (see Figure|3J). There are plain vertices, anti-vertices, and empty vertices, 
as well as plain vertical edges (vedges), anti-vedges, and plain horizontal edges (hedges). For 
now, plain components are given weight 1, anti- vertices and anti-vedges are given weight -1, and 
empty vertices have weight zero.^ 

Plain Anti Empty 
Vertices • O □ 

Vedges 

Hedges 



Figure 3. Basic components of a signed graph. 



Define a signed matching of a graph G to be a collection of non-adjacent edges of G and all 
the vertices not incident with those edges. The weight of the matching is the product of the 
weights of the components in the matching. Note that this implies that a matching that involves 
an edge one or both of whose endpoints is an empty vertex has vanishing weight. Matchings of 
weight zero "do not count" for purposes of weighted enumeration, so in practice it helps to think 
of empty vertices as being needy. For any matching of nonzero weight, empty vertices need to 
have one of their incident edges included. 



□- 



-O- 



□- 



-o- 



FlGURE 4. Pictures of graphs G(2,2), G(3,-3), and G*(2,2). 

Define a regular signed graph, G(m, n) = Pm x Pn iu the following way. If n > 0, G{m,n) 
is just a rectangular grid graph, i.e., a graph with mn plain vertices arranged into m rows of n 
vertices each, with plain edges adjoining horizontal and vertical neighbors. Define the conjugate 
graph of G, denoted G*, to be the graph obtained by replacing all of the vertices and vedges 
of G with their anti-counterparts (assuming that an anti- anti- vertex is a plain vertex, etc.) but 
leaving the hedges alone. Under this definition, (G(m,n))* = G*(m,n) is a graph with mn 
anti- vertices arranged into m rows of n vertices with hedges adjoining horizontal neighbors and 
anti-vedges adjoining vertical neighbors. Then for n < 0, G{m,n) is defined to be a copy of 
G*{m,n) with a column of m empty vertices on the left and m empty vertices on the right. 
Each empty vertex is connected (with a hedge) only to the anti-vertex horizontally adjacent to 
it (see Figure 

For fixed to > and any integer n, we define the number M(to, n) — M[G{m^n)) to be the 
sum of the weights of the signed matchings of G{m,n). In general, M{G) is the sum of the 
weights of all of the signed matchings of G. Some examples are given in Figure [S] 

Having made all our definitions in terms of rectangular grid graphs, we would now like to 
extend them to apply to generalized rectangles. For any finite graph G, define the generalized 
rectangle graph G x P„ with n > as follows: Picture n copies of G lined up side by side. Each 

^Later we will attach formal variables to edges and vertices, but empty vertices will still have weight zero. 
Thinking of monomer-dimer coverings, empty vertices represent positions in the lattice that may not be occupied 
by a monomer. 
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G(2,2) 



+1 



+1 



+1 



+1 



+1 



+1 



+1 



□ o o o- 



n o 



n o o o- 

□ o o o 



□ o o o- 



□- 



□- 

o 



□- 



-o 



-o 



G(3,-3) 



+1 



+1 



+1 

M(G(3,-3)) 



+7 

M(G(2,2)) 



Figure 5. Pictures of nonzero signed matchings of G(2, 2) and G(3, —3). 



vertex in one copy of G is connected with an edge to its image in the copies of G to the left and 
right of it. Define all the edges within a copy of G to be vedges, and the edges between copies of 
G to be hedges. With this idea, we can define G x P_„ to be G* x P„ with m {— the number of 
vertices in G) empty vertices to the left and to the right of it. Each empty vertex is connected 
to exactly one of the vertices of G* with a hedge. We shall see that these definitions are the 
right ones that will enable us to prove statements of reciprocity for more than just grid graphs. 

3. Adjunction 

Although it seems there are separate definitions for G(m,n) (resp. G x P„) when n > and 
when n < 0, we will show they are actually the same object. In Propp's words, they fit together 
"seamlessly." If they fit together this way, then it follows that for any G, the sequence M{G x P„) 
satisfies a linear recurrence of finite degree with constant coefficients. We will omit the proof as 
it is given in . The method of proof utilizes transfer matrix methods and the Cayley-Hamilton 
theorem, and it requires only that the family of graphs satisfy the property detailed below. 

We begin by making a natural observation about the half-infinite sequence M(G(m, n)), ti > 0, 
and then we will prove that it actually holds for any n. In doing so, we will be able to establish 
that we do indeed have an appropriate definition for G{m,n) (resp. G x Pn) with n < 0, and 
all our desires for the doubly infinite sequence — integrality, uniqueness, and reciprocity — will be 
fulfilled. 



A RECIPROCITY THEOREM FOR MONOMER-DIMER COVERINGS 



7 



G(2,2) 



G(2,3) 




G(3,-3) 

□ o o o- 



□- 



-o- 



-o- 



□ o o o □ 



□ o o o- 



n o o 



-o o o □- 

7 



G(3,2) 



G(2, 2)G(2, 3) = G(2, 5) G(3, -3)G{3, 2) ~ G(3, -1) 

Figure 6. Pictures of two pairs of adjoined graphs. 



We first notice that for positive ni,n2, we can draw the graph G{m,ni + 712) by placing 
G{m, ni) and G^m, side by side and using hedges to connect the rightmost m vertices of one 
to the leftmost m vertices of the other, as shown in Figure |S| Call this operation adjunction. 
Generally, define the adjunction of two graphs H = G x P„j , H' — G x Pn^ to be a new graph 
formed by connecting H to H' with hedges, matching up corresponding vertices. In particular: 

• If JT-i, 77.2 are positive, connect every vertex in the rightmost copy of G in H to its image 
in the leftmost copy of G in H' using a hedge. 

• If 7T.1 is positive, 77,2 negative, join every vertex v in the rightmost copy of G in to 
an empty vertex on the left side of H' so that v is connected to v* (the anti-vertex 
corresponding to v) in the leftmost copy of G* in H' by the path hedge-empty vertex- 
hedge. (Similarly if rii is negative and 7t2 is positive.) 

• If 771, 772 are both negative, join every vertex in the rightmost copy of G* in H to its 
image in the leftmost copy of G* in H' by connecting their adjacent empty vertices with 
a hedge. 

We write the adjunction of G x P„^ and G x P„2 as (G x P„J(G x Pnj). 

Returning to case of rectangular grid graphs where 771,772 are positive, G(777, 7ii)G(777, 772) = 
G(777 , 77i -t-772), and naturally, M(G(777 , 77i)G(777, 712)) = M{G(m,ni + 772)). This is true more 
generally, as stated in the following: 

Theorem 1 (Adjunction). 

M((G X P„J(G X P„J • • • (G X P„J) ^ M(G x P„, 

-1-712-1 hrifc ) 

for all integers 7ti, . . . , 77^. In particular, 

M{G{m, ni)G{m, 772) • • • G{m, 77^)) = M{G{m, ?7i -|- ?72 + ■ • • + f^fe)) 

for all integers 77i, . . . , 77^. 

The proof of Theorem ^ will require the following two lemmas. 

Lemma 1. Let G be a signed graph where two vertices, a and b, are connected by the path 
hedge-empty vertex-hedge-empty vertex-hedge. Then G has same number of signed matchings as 
the graph G' , where G' is identical to G except that a and b are connected with one hedge. 
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Figure 7. Proof of bijcction between G and G' in Lemma 1. 

Proof: The claim made by Lemma 1 is most easily seen in Figure|Z|where we see the immediate 
bijection. When a matching of G does not contain the hedge joining the two empty vertices, 
then the matching must contain the two edges that join these two empty vertices to a and b. 
This corresponds to a matching of G' that contains the hedge joining a and b. 

On the other hand, any matching of G where the hedge between the empty vertices is present 
is a matching in which a and b are connected outwards. Such a matching corresponds to a 
matching of G' that does not contain the hedge joining a and b. □ 

Lemma 2. Let G be a finite graph with all plain components. Then let H be the graph (G x 
Pni ) (G X ) for ni > > 712 . Then H has the same number of signed matchings as the graph 
H'={G X P„,_i)(G X P„,+i). 

Proof: We begin by examining what happens at each vertex; see Figure |S1 Consider the plain 
vertex a of G connected to the anti-vertex a* of G* by the path hedge-empty vertex-hedge. 
There are essentially two types of matchings of H . We will say type 1 matchings include the 
hedge connecting a to the empty vertex and type 2 matchings include the hedge connecting a* 
to the empty vertex. Most of the matchings of type 1 will cancel with most of the matchings of 
type 2. 

Given a matching of type 1, there are only three cases for what can happen the anti- vertex 
a* of G*: it can be isolated, it can have an anti-vedge incident with it, or it can have a hedge 
incident with it. Likewise for any matching of type 2, the vertex a can be isolated, incident with 
a vedge or incident with a hedge. If there are k signed matchings of type 1 where a* is isolated, 
then there are —k signed matchings of type 2 where a is isolated. Similar cancellation occurs 
between matchings of type 1 where a* has an anti-vedge and matchings of type 2 where a has a 
vedge (corresponding to the anti-vedge incident with a*). 

The only remaining cases are those of type 1 and type 2 where a and a* both have hedges. 
We claim that these matchings are in bijection with the graph where a, a*, all their incident 
vedges, and the hedges and empty vertex between a, a* are all replaced with one empty vertex. 
The correspondence is shown in Figure |H1 Since a was any vertex of G, this can be done for 
every vertex of G and the lemma holds. □ 

With the lemmas proved. Theorem is not difficult to show. 

Proof of Theorem 1: It is easily verified that M((G x Pq){G x P„)) = M(G x P„), so we 
may assume all the Ui are non-zero. We have that M{{G x Pn^){G x Pna)) — M{G x Pn^+^j) 
whenever ni and 712 are negative by applying Lemma 1 to each place where adjoining takes 
place. The case ni,7i2 > is even easier. When ni and 712 are of opposite sign, we can apply 
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H H' 

G G* 



+/- 




Type 2 



Figure 8. Proof of Lemma 2. 

Lemma 2 repeatedly. Upon each application of Lemma 2 we change neither the number of signed 
matchings nor the difference between the number of vertices and the number of anti- vertices. □ 

4. Reciprocity 

4.1. Combinatorial Reciprocity. The combinatorial statement of reciprocity is rather obvi- 
ous given the definition of G x for n < and the neediness of empty vertices: 

Theorem 2 (Reciprocity (I)). 

M(G X F_„_2) = M{G* X P„). 

In particular, 

M(G(m,-n-2)) = M(G*(m, n)). 

Proof: On the left side of FigurelTUlwe have G(m, —n — 2). But since all the outside edges are 
forced and have weight 1, the number of signed matchings of G(m, —n — 2) is clearly equal to 
the number of matchings of G*{m,n). The situation for generalized rectangles is identical. □ 

This statement of reciprocity is the natural extension of the property of dimer coverings 
observed by Propp (0)- Indeed, if we were to consider only the perfect matchings of grid 
graphs, then the statement would be the same: 



7V(G(m, -71 - 2)) = iV(G*(TO,7i)). 
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G{in, n2) 



Lemma 2 

ni > > n2 




Figure 9. Picture of Theorem using lemmas. 
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Figure 10. Picture of reciprocity for grid graphs. 
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G*(m,n) 
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It is known that any dimer covering can be obtained from any other dimer covering by local 
moves that leave the parity of vedges unchanged. That is, the sign of every perfect matching 
of G*{m,n) is the same, so \N{m,n)\ — \N*{m,n)\. With less-than-perfect matchings, this is 
clearly not the case. 

Perhaps the strength of this result is more obvious when we give all the edges and vertices 
formal weights and derive a statement of reciprocity in terms of a generating function. Though 
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we could do so for any finite graph of tlic form G x P„, we will only derive an explicit formula 
for G(m, n). 

4.2. The Generating Function. For convenience, let the weights of the vertices of G{m,n) 
for TT, > be indexed by N x N, read from left to right and top to bottom. Notice that the top 
left vertex is 2:1^1 and the bottom right vertex is Zm,n- Then Xi^ is the weight given to the hedge 
between Zij and Zi,j+i, and ytj is the weight given to the vedge between Zij and Zi+ij- To get 
the weights for vertices of G{m, —n), we begin by giving the top right anti- vertex weight 2:1^0 
and we proceed analogously. The weights of G(2, 3) and G(2, —5) are shown in Figure [TTl 




□ □ 

yi,-4, I yiys 1/1,-2 I vl,-! yi.o 

□ — — — — — □ 



Figure 11. Picture of graphs weighted with formal variables. 

The weight of a matching with formal variables is not exactly like the weight as defined before, 
as is seen in the following: 

Definition 1. If ^ is a matching of G(rn,n), define the weight of ii, as follows (where 

sgn(t) = 1 if t is positive, —1 otherwise): 



win) = 



n 



X. 



hedges in /i vcdgcs in fj. 



vertices in /i 



n 



X. 



all hedges in G(m,n) 



(l-sgn(i))/2 



Notice that for n < 0, there is a monomial denominator in the weights of matchings. By our 
definition, the matchings shown in Figure [T^ have weights 



a;i, 22/1, 12^2, 2^2, 3 and 



-12:^1, -2a;i, -42:2, -ia;2,-2a;2, -32^2, -4 



Figure 12. Picture of matchings weighted with formal variables. 
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With this new definition of weights, we define the matching polynomial. In the literature, 
there are conventions for describing the matching polynomial of a graph. Our definition does 
not adhere strictly to these conventions, though it is similar to the partition function for the 
monomer-dimer model, with distinct formal variables for all monomers and dimers (see 

Definition 2. The matching polynomial of G{m,n) , fn{xij,yij, Zij), is given by 

fn{Xij,yiJ,Z^j) = ^ w{fi) 

all matchings of G(m,n) 

where i = 1, . . . ,m, j = 1, . . . ,n if n > 0, j = 0, —1, . . . ,n if n < 0. 

Strictly speaking, the matching polynomial is a Laurent polynomial: a polynomial in the 
variables Xij, x~j, yij, Zij. Notice if we set all of the weights equal to 1 we get /«(!, 1, 1) — 
M{G{m, n)), the number of signed matchings. Similarly if we set all the Xij,yij equal to 1 but 
set all the Zij equal to 0, then /n(l, 1, 0) — N{G{m, n)), the number of signed perfect matchings. 

By construction the polynomials /„ will satisfy a linear recurrence very similar to that which 
governs M{G{m, n)). For example, if to = 1, then we get /„ = zi^nfn-i + xi^n~ifn-2 for all n. 
Likewise, there is a link between the polynomials /„ and f-n-2 for n > 0, seen most easily if we 
let Xij = X, yij — y, Zij — z for all j, j: 

(1) /n(a:, -2/, ~z) = x™("+i)/_„_2(a:, z). 

As mentioned in Section 3, for any fixed to, the numbers M{m,n), n G Z satisfy a linear 
recurrence of finite degree with constant coefficients. Therefore there is a rational generating 
function for the number of weighted matchings, 

F„i{t, X, y,z) = fn{x, y, z)t". 

n>0 

Given (1) it is an easy exercise (see for example 0, ch. 4) to state a reciprocity theorem for the 
generating function: 

Theorem 3 (Reciprocity(II)). 

x™-t^Fm{t, X, y, z) = -F,n (j^'X' -V' ■ 
This completes the main goal of the paper. 



5. More on Linear Recurrences and Reciprocity 

Though we have only proved theorems here for graphs of the type G x P„ (cylinders of fixed 
circumference and varying height for example), we have also been able to use a modified form 
of adjunction to build a model for the honeycomb graph of fixed height and width n G Z. Other 
graphs that we examined but were unable to apply our methods to include cylinders with fixed 
height and varying circumference, Mobius strips, tori, and projective planes. We do however feel 
that there may be a way of dealing with such graphs. 

In his unpublished paper [Sj, David Speyer developed a matrix method for encoding perfect 
matchings of a graph. With this method he was able to state theorems about recurrences and 
reciprocity for a broader range of graphs than those we have handled here. In particular, he was 
able to make a statement about the Mobius strip and projective plane. An adaptation of his 
method to general matchings seems promising, though experimentation has shown it to be less 
than straightforward. 



A RECIPROCITY THEOREM FOR MONOMER-DIMER COVERINGS 



13 



6. Acknowledgements 

We'd like to thank Jim Propp for introducing us to this problem and for getting us started 
with the early investigations. Without his REACH program to bring the authors together, this 

work could not have been done. Thanks also to Harvard University, Brandeis University, the 
University of Massachusetts-Boston, the University of Wisconsin, the NSF, and the NSA. 

References 

[1] Y.B. Choe, J.G. Oxley, A.D. Sokal, D.G. Wagner, Homogeneous multivariate polynomials with the half-plane 
property, preprint. 

[2] O.J. Heilmann, E.H. Lieb, Theory of monomer-dimer systems. Comm. Math. Physics, 25 (1972), 190-232. 
[3] P.W. Kasteleyn, The statistics of dimers on a lattice, Physica 27 (1961), 1209-1225. 

[4] C. Kenyon, D. Randall, and A. Sinclair, Approximating the number of monomer-dimer coverings of a lattice. 

Journal of Statistical Physics 83 (1996), 637—659. 
[5] J. Propp, A reciprocity theorem for domino tilings, The Electronic Journal of Combinatorics 8 (2001), R18. 
[6] D. Speyer, Unpublished note on transfer matrices. 

[7] R. Stanley, Enumerative Combinatorics, Volume I, Cambridge University Press, 1997. 
[8] R. Stanley, On dimer coverings of rectangles of fixed width. Disc. Appl. Math. 12 (1985), 81-87. 
[9] H.N.V. Temperley and M.E. Fisher, Dimer problem in statistical mechanics - An exact result. Philosophical 
Magazine 6 (1961), 1061-1063. 



